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ABSTRACT 

We show a possible way to measure the column density of free electrons along the light path, the 
so-called Dispersion Measure (DM), from the early [~ 415(^/1 GHz) _2 (DM/10 5 pc cm -3 ) s] radio 
afterglows of the gamma-ray bursts. We find that the proposed Square Kilometer Array can detect 
bright radio afterglows around the time ~ 10 3 (^/160 MHz) -2 s to measure the intergalactic DM 
(> 6000 pc cm -3 at redshift z > 6) up to z ~ 30, from which we can determine the reionization 
history of the universe and identify the missing warm-hot baryons if many DMs can be measured. 
At low z, DM in the host galaxy may reach ~ 10 5 pc cm -3 depending on the burst environment, 
which may be probed by the current detectors. Free-free absorption and diffractive scattering may 
also affect the radio emission in a high density. 

Subject headings: gamma rays: bursts — intergalactic medium — radio continuum: ISM 



1. INTRODUCTION AND SUMMARY 

Many observations support a massive star origin for 
long-duration Gamma- Ray Bursts (GRBs). Although 
the nature of the central engine that accelerates the GRB 
jets remains unknown, the afterglows are successfully fit- 
ted by the synchrotron shock model for probing the jet 
properties and burst environment (e.g., Meszaros 2002). 

On the other hand, GRBs may be useful for probing 
high redshift z (e.g., Miralda-Escude 1998; Barkana & 
Loeb 2003; Inoue, Yamazaki, & Nakamura 2003). Their 
high luminosities make them detectable even at z ~ 100 
(Lamb & Reichart 2000), while the X-ray and infrared 
afterglows are observable up to z ~ 30 (Ciardi & Loeb 
2000; Gou et al. 2003) and provide their z through the 
Lya break (Lamb & Reichart 2000) or Fe lines (Meszaros 
& Rees 2003), in contrast to galaxies or quasars which 
are likely dimmer at higher z. About 25% of all GRBs 
detected by the upcoming Swift satellite are expected to 
be at z > 5 (Bromm & Loeb 2002). High z GRBs may 
have already been detected by BATSE based on the em- 
pirical luminosity indicators (Fenimore & Ramirez-Ruiz 
2000; Norris, Marani, & Bonnell 2000; Murakami et al. 
2003; Yonetoku et al. 2003), for which some theoretical 
explanations exist (e.g., Ioka & Nakamura 2001a). The 
first generation stars could be very massive (Abel, Bryan, 
& Norman 2002; Bromm, Coppi, & Larson 2002; Omukai 
& Palla 2003), so that they may end as brighter GRBs. 

In this Letter, we show that free electrons along the 
light path cause distortion in the spectrum and light 
curve of an early radio afterglow, from which we can 
measure the column density of the free electrons, the so- 
called Dispersion Measure (DM). 2 This is because, in a 
plasma with an electron density n e , an electromagnetic 
wave with frequency v(~^> i/ p ) is delayed relative to in a 
vacuum by a time 

f dl vl t v \ - 2 / DM \ , s 
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2 Inoue (2003) has also independently carried out a similar study. 



where v v = (n e e 2 /irrrie) 1 / 2 = 8.98 x 10 3 nl^ 2 Hz is the 
plasma frequency and DM = J n e dl (Rybicki & Light- 
man 1979). DMs of pulsars are well studied by using the 
arrival time of pulses at two or more frequencies (Tay- 
lor, Manchester, & Lyne 1993). Our method to measure 
DMs of GRBs is somewhat different since we do not as- 
sume simultaneously emitted pulses but the afterglow 
model. We can determine DM only from a single-band 
light curve around the time t ~ At (see The former 
method was suggested by Ginzburg (1973) and Palmer 
(1993) before the discovery of the afterglows. 

At z > 6, DMigp due to the intergalactic plasma 
(IGP) will be > 6000 pc cm" 3 [correspondingly At > 
10 3 (i//160 MHz)" 2 s], and probably dominate DM G due 
to the Galactic plasma and DMhost due to the plasma 
in the host galaxy (see 0- DMtqp as a function of z 
varies depending on the reionization history of the uni- 
verse since recombined electrons provide no DM (see 
Thus, from DMs of GRBs at various z and directions, 
we can determine the reionization history, and possibly 
even map the topology of the ionized bubbles. The reion- 
ization history is now actively investigated but remains 
unclear (Miralda-Escude 2003). The analysis of the Lya 
spectra in the highest z quasars suggests that the reion- 
ization ends at z ~ 6 (Fan et al. 2002), while the WMAP 
polarization data imply a much higher reionization red- 
shift z ~ 17 ±5 (Kogut et al. 2003; Spergel et al. 2003). 
In 21 we show that the proposed Square Kilometer Ar- 
ray (SKA) can detect bright afterglows at t ~ At to 
measure DMtqp up to z ~ 30. In addition, observations 
of DMiqp could also identify the 'missing' baryons that 
have not been detected at z < 1 (see 

At low z, DMhost in the host galaxy may reach ~ 
10 5 pc cm -3 and dominate DMq and DMtqp depend- 
ing on the as-yet-unknown burst environment (see SjSJl. 
Since we may target low z GRBs, even current detec- 
tors could detect afterglows at t ~ At to constrain 
DMi los t and the GRB environment. Throughout we 
adopt a ACDM cosmology with (Q m ,Q\,ili,,h,as) — 
(0.27,0.73,0.044,0.71,0.84) (Spergel et al. 2003). 

2. EXPECTED DISPERSION MEASURE (DM) 
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Fig. 1. — The Dispersion Measure (DM) due to the intergalactic 
plasma, DMigp, is plotted as a function of redshift z by the solid 
line. The dispersion of DMigp due to density fluctuations is shown 
by the shaded region. 0.5 X DM™ a *(z) and 0.1 X DM™ a *(z) are also 
shown by the dotted lines where DM™ a *(z) = DM™ a */(l+z) is the 
maximum DM due to the plasma in the host galaxy in equation 
The dashed lines are DMigp for a sudden reionization at 
^rcion = 6, z rc ion = 17 and the case when the reionization occurred 
twice at z rc i on = 6, 20 with a full recombination at z = 15. 



Galactic DMq: The distribution of free electrons in our 
Galaxy is relatively well known from pulsar DMs (Tay- 
lor & Cordes 1993). DMq has its maximum (minimum) 
in the direction parallel (perpendicular) to the Galactic 
plane, ranging from 



DMg m 



30 pc cm 



to DMg ax 



10 3 pc cm" 3 (Taylor et al. 1993; Nordgren, 
Cordes, & Terzian 1992). 

Intergalactic DMigp'- We first estimate DMiqp as- 
suming that all baryons are fully ionized and homo- 
geneously distributed, so that the free electron den- 
sity evolves as n e = (3H 2 rib/8irGm p ) (1 + z) 3 . 
Then the observed dispersion delay at an observed 
frequency v for a source at z is given by At = 
f* dz\dt/dz\(l + z)v 2 /2[v(l + z)} 2 = 415(^/1 GHz) -2 



1/2 



(DMigp/10 5 pc cm -3 ) s, where \dt/dz\ 

z)h [n m (i + zf + n A ] 

3cHo^tb 



(1 



and 



DMigp = 



87tGto 



(1 + z)dz 



p Jo [fi m (i + z) 3 + n A ] 



1/2' 



(2) 



From FigureH DM IGP >DMg in at z > 0.03. 

More than half of all baryons have not been detected at 
z < 1 (Fukugita, Hogan, & Peebles 1998). Most of such 
'missing' baryons may reside in a warm-hot IGP, which 
is difficult to observe because of its high temperature 
(10 5 -10 7 K) and low density (moderate overdensity 10- 
40) (Cen & Ostriker 1999; Dave et al. 2001). If so, an 
appreciable fraction of DMiqp Qt 500 pc cm" 3 ) comes 
from the missing baryons at < z < 1 (see Figure QJ. 
Thus observations of DMiqp could be a direct detection 
of the missing baryons. 

DMhost in host galaxy: A host galaxy has been found 
in most GRBs (Bloom, Kulkarni, & Djorgovski 2002). 
If the host galaxies are similar to ours, they have 



DMi los t ~ DMq. At high z, the maximum and minimum 
DM host may evolve as ~DMg ax (M/10 12 M Q )(l + z) 3 and 
~DMg in (A'//10 12 M Q ) 1 / 3 (l + zf, respectively, for a host 
galaxy with the halo mass M based on the hierarchi- 
cal galaxy formation (Ciardi & Loeb 2000). Even at 
low z, DMhost may be high if GRBs arise from the 
star forming regions, as expected in GRBs resulting 
from massive stellar collapses. The hydrogen column 
density is as high as 10 22 -10 23 cm" 2 in giant molec- 
ular clouds (Galama & Wijers 2001), and the prompt 
and afterglow emission with energy E lon erg can ion- 
ize the ambient density rihost cm" 3 out to a distance 
of d ~ ^0n~^ t 3 E^l 52 pc (Perna & Loeb 1998), where 
the convention Q = 10 X Q X is used. Thus all the ambi- 
ent hydrogen may be ionized depending on its density. 
In this case, DMhost ~ 10 3 -10 5 pc cm" 3 , which may be 
time-dependent (Perna & Loeb 1998). A high DM is also 
expected in the 'supranova' model (Vietri & Stella 1998), 
in which a supernova (SN) occurs weeks before the GRB. 
A SN remnant shell with mass M g and velocity v cm 
s _1 reaches at R ~ 10 16 vgi7 cm in t s, which provides 



1.5 x lOSg 2 if 2 M 33 pc cm" 3 



DMhost ~ M/^R 2 m p 

Anyway, the maximum DMhost allowed by the present 
observations is about the inverse of the Thomson cross 
section <jt, 



r>M max 



4.87 x 10 5 pc cm" 



(3) 



for GRBs and afterglows to be observed without scat- 
tering. Correspondingly, Aihost ~ 2022(^h ost /l GHz) -2 
s from equation 0). At high z, the observed disper- 
sion delay and frequency are At = (1 + z)Athost and 
v = ^host/(l + z), respectively. Thus it is convenient to 
use DMhost(^) = DMhost/(l+2 ; ) to compare it with DMq 
and DMigp. From Figure Q] DMhost (z) declines with in- 
creasing z while DMigp grows. Therefore, GRBs at high 
(low) z are suitable for measuring DMiqp (DMhost)- 

3. DM FROM DISPERSED AFTERGLOWS 

We now study the effects of the DM on the afterglow 
emission. For illustrative purposes we fix DM demanding 
the observed dispersion delay of At = 10 3 (i^/l GHz)" 2 s 
in this section. We consider the most simple but standard 
afterglow model (Sari, Piran, & Narayan 1998), which 
reproduces the observations very well. 

The afterglow spectrum for the forward shock is well 
approximated by four power law segments with breaks 
at cooling frequency z/ c , typical frequency v m and self- 
absorption frequency v a . In the fast cooling case, [y a < 



)u c < v n 



the spectrum is F v cx v 



.1/3 ,,-1/2 



-p/2 



from low to high frequencies, and has a maximum flux 
-fl/,max at v = v c . In the slow cooling, (y a <)v m < v c , it 
is F v cx v 2 , h> l / z , i/"( p_1 )/ 2 , v~ v l 2 , and peaks at v = v m . 
Assuming an adiabatic shock, we have the observed break 
frequencies and maximum flux as 



v c = 2.h x 10 ld e 



13 -3/2 p-l/2 -1,-1/2 



B E 52 



1-1/2 



Hz, 



u m =4.2 x 10™e)/ 2 e 2 e g 2 Ell 2 t^ 2 (l + zf' 2 Hz, 

i^.rnax = 1.1 X 10 5 4 /2 £ 52 n 1/2 £28 2 (l + z) ^Y, 

where E erg is the isotropic equivalent shock energy, n 
cm" 3 is the constant surrounding density, e e (es) is the 
electron (magnetic) energy fraction, t s is the observer 
time, D cm is the luminosity distance and g = (p — 



Dispersion Measure from GRB Afterglows 
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Fig. 2. — The spectra (upper panel) and light curves (lower 
panel) of afterglows are shown with and without the dispersion de- 
lay At = 10 3 (u/l GHz) -2 s by the solid and dashed lines, respec- 
tively. We adopt E = 10 52 erg, n = 1 cm -3 , e e = 0.1, e B = 0.01, 
p = 2.2 and z=0.5. For these parameters, the afterglows in this 
figure are in the slow cooling regime. 



2)/(p— 1). At v < v a , synchrotron self-absorption limits 
the flux below the blackbody emission with the electron 
temperature (e.g., Sari & Piran 1999; Kobayashi & Sari 
2000). This is given by 

F„, BB = 27w 2 77e m e (R±/D) 2 (1 + z), 

where v is the observed frequency, 7 is the Lorentz factor 
of the shocked fluid, R± = 4 7 ct is the observed size of 
the afterglow and 7e is the typical Lorentz factor of the 
electrons emitting at v. A reverse shock is also formed, 
but its temperature is lower because of its higher density, 
making its emission dimmer at v < v a . 

Once we calculate the spectrum as a function of time 
F v (i), we can obtain a dispersed afterglow as F v [t — 
At{v)] where At(v) is the dispersion delay. From Fig- 
ure |3 the spectrum has a cutoff at v ~ l(t/AtGHz) -1 ^ 2 
GHz when At = A£ G HzO/1 GHz)" 2 where we are us- 
ing AtcHz = 10 3 s. The light curve deviates from the 
power law at t <~ At. Measurements of these features 
yield At and hence DM with equation Q. In practice, 
DM may be treated as an additional fitting parameter. 
Only a single-band light curve will suffice to measure 
DM, but a multiband observation will be useful to ex- 
clude other interpretations, such as variable light curves 
seen in GRB 021004 and GRB 030329 (e.g., Uemura et 
al. 2003) and microlensing (Loeb & Perna 1998; Ioka 
& Nakamura 2001b), by using the frequency dependence 
At tx v~ 2 . 

4. PROSPECTS FOR PROBING THE REIONIZATION 

DMigp depends on the reionization history (see Figure 
nj. DMigp is a constant ~ 6000 pc cm~ 3 at z > z rcion 
if a sudden reionization occurred at z ro i on = 6, while 
DMigp ~ 12000 pc cm -3 if z rcion = 17. We also show 
the case when the reionization occurred twice at z rc i on = 
6, 20 with a full recombination at z = 15 (Cen 2003a, b; 
Wyithe & Loeb 2003). Contrary to the Lya absorption 



Fig. 3. — The afterglow flux at a frequency v = 160 MHz and 
an observer time t = 10 3 s is shown as a function of redshift z 
by the solid lines. We adopt e e = 0.1, eg = 0.01 and p = 2.2 
with E = 10 52 , 10 54 , 10 52 erg and n = 1, 10" 2 , (1 + z) 4 
cm -3 for typical, bright, dim afterglows, respectively. The rms 
sensitivity of the Square Kilometer Array for an integration 
time 10 3 s and a band width Au ~ 0.5u is plotted by the 
dashed line. The expected flux of the host galaxy F^ OBt ~ 
25[(1-M>/1 GHz]" ' 75 (SFR/1 M yr" 1 ) (D/l Mpc) _2 (l + 
z) Jy is shown by the dotted lines where SFR is the star formation 
rate (Yun & Carilli 2002). 



(Inoue et al. 2003), DM can determine the reionization 
history even if the neutral fraction is > 10~ 5 in the first 
reionized era. If DMhost(^) >DMiqp, the mean DM of 
the GRBs will be above the line of DMigp m Figure 

Can we measure DMs of high z GRBs ? We 
should start a radio follow-up within ~ At ~ 
972 (j//160 MHz)" 2 (DM/6000 pc cm -3 ) s, since DM 
distorts the light curve at t ~ At (see Figure EJ- 
This may be possible since Swift will send a 1-4 ar- 
cmin GRB position to the ground within ~ 1 min (see 
http://swift.gsfc.nasa.gov/). The dispersion delay At is 
longer at lower frequencies, but the flux decays as F v oc 
v 2 at v < v a - Here we adopt v — 160 MHz assuming a 
response time of t ~ 10 3 s. In Figure|3 the afterglow flux 
at v = 160 MHz and t = 10 3 s is plotted as a function 
of z. For the typical parameters, E — 10 52 erg, n = 1 
cm~ 3 , e e = 0.1, €b = 0.01 and p = 2.2, the flux is too dim 
to be detected even by SKA, whose rms sensitivity is 

0.3 (Av/B0 MHz)" , , 
0.5v, an integration time r 



^ 2 t 3 X ^ 2 /iJy for a band width Av 



10 3 s and A eS /T sys = 2 x 
1 (see http://www.astron.nl/skai/science/). 
However some GRBs have larger energy E ~ 10°* erg 
and lower density n ~ 10 -2 cm~ 3 (Panaitescu & Ku- 
mar 2002). From Figure E3 such GRBs can be detected 
as ~ 7cr events even at z = 15 by SKA, and dominate 
the host galaxy emission even if SFR ~ 100A/ Q yr _1 at 
z >8. Note that the density around the first stars could 
be 10~ 2 < n < 1 cm~ 3 because of strong radiation pres- 
sure from the central massive star (Gou et al. 2003). On 
the contrary, the density could evolve as n oc (1 + z) 4 for 
a fixed host galaxy mass based on the hierarchical galaxy 
formation (Ciardi & Loeb 2000). If n - (1 + z) 4 cm~ 3 , 
we will never detect radio afterglows at high z. 
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5. DISCUSSIONS 

Free-free absorption (FF) may be important when 
the ambient density nh os t cm -3 is high. The 
host galaxy is optically thick to FF at v < 
2(DM host /10 5 pc cm- 3 ) 1 /2T 4 ' 3/4 n^ a(1 + z yi qjj z 

where T is the plasma temperature in K (Rybicki & 
Lightman 1979). In this case we cannot measure DM 
but can probe the burst environment. We may neglect 
FF in IGP at v > 1.6/i 3 / 2 rj 6 ^ 2 ^ m ^ 4 iT 4 _3/4 (l + z) 5 / 4 kHz 
where z > 1 (Rees 1978). 

Large-scale density fluctuations, as well as the accu- 
mulation of gas in discrete galaxies, cause a dispersion 
in DMigp- For example, if the average number of in- 
tervening gas clumps with a DM of £ is TV, the dis- 
persion ((ADM) 2 ) is ~ NT, 2 if we assume a Poisson 
statistic. We may roughly estimate ((ADM) 2 ) using the 
Press-Schechter theory. Kitayama & Suto (1996) de- 
rived the comoving number density of halos that form 
with mass M ~ M + dM at time Zf ~ Zf + dzf and 
are observed at z, F(M, zf, z)dMdzf. If the halo mass 
is M < 10 12 M©, the inside gas can cool to have a ra- 
dius r d ~ 0.035r vir (e.g., Ciardi & Loeb 2000), oth- 
erwise rd ~ r v ; r , where r V j r (M, zf) is the virial ra- 
dius of the halo. Then, DM of one halo is about 
T,(M,Zf,z) — MQb/2'Kr 2 .m p Q m (l + z), and the disper- 
sion of DM due to halos in a logarithmic mass interval 
is about ([ADM(M)] 2 ) - / cdtj dz f Trr 2 MF(l + zfY?. 
We calculate max([ADM(M)] 2 ) and show it as the 
shaded region in Figure ^ It takes the maximum when 
M ~ 10 12 M Q and is about ([ADM(10 12 M Q )] 2 )/DM 2 ~ 
(68%) 2 at z = 1, though the mean number of interven- 
ing halos is ~ 10 -3 . Unvirialized objects may have less 
contributions since ([ADM(2 x 10 14 M Q )] 2 )/DM 2 - 8/t" 1 



Mpc/6.6 Gpc ~ (4%) 2 at z = 1 for density fluctuations 
of radius 8ft. -1 Mpc. 

At v < v s — 10 GHz, diffractive and refractive scin- 
tillation due to the Galactic plasma take place (Good- 
man 1997; Walker 1998). Although the modulation 
index (rms fractional flux variation) is unity m s = 1 
for diffractive scintillation, the decorrelation bandwidth 

is very narrow /S.v s jv — {v/v s ) 17 / 5 ~ 4 x lO -4 ^ 7 ^ 5 
at low frequencies. Then m s will be reduced to ~ 
yj Ai/ s /u ~ 2^g 7 ^ 10 % for a broadband observation. 
It will be also reduced when the integration time is 
longer than the scintillation one ~ 700i^g /5 s. For re- 
fractive scintillation, m s = (v/v s ) 17 ^ 30 ~ 27vg 7 ^ 30 %. 
Thus Galactic scintillation may be quenched at v < 1 
GHz with Af ~ v. On the other hand, scattering 
in the host galaxy may cause temporal broadening of 
the early radio afterglows by a time ~ l/2irAv a ~ 

10^ 22/5 (DM host /10 5 pc cm-Y/Bn^dxCl + z)" 17 / 5 
s, if the scattering measure scales as cx In 2 (Goodman 
1997), where d pc is the size of the ionized region. 

Some GRBs may have the wind environment (Cheva- 
lier & Li 2000). Inverse Compton may prolong the fast 
cooling regime (Sari & Esin 2001), in which an additional 
spectral segment F v oc v 11 / 8 may appear (Granot, Piran, 
& Sari 2000). In addition, photons accumulating around 
the shock front may make the cooling 'ultrafast' (Ioka 
2003). These are interesting future problems. 

We thank F. Takahara, J. Yokoyama, K. Omukai and 
M. Matsumiya for useful discussions. This work was sup- 
ported in part by the Monbukagaku-sho Grant-in-Aid 
No. 00660. 
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